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I. INTRODUCTION 


Satisfactory scenarios of orbit and spectrum allocations for 
communications satellites must achieve sufficient protection for the 
desired signals relative to interfering ones. Since the protection 
ratio is a highly non-linear function of the satellite locations and the 
frequencies to be assigned, it is natural to formulate the optimization 
of these assignments in terms of non-linear programming techniques 
[1-4]. The resulting codes place heavy demands on computational 
resources [5]. Here we shall introduce an intermediate step: from an 

estimate of the required single-entry protection we shall first 
calculate the required satellite separations and then use these as 
constraints in a linear-programming optimization, which should require 
less computational effort. To do so, it is first necessary to derive 
quantitative relationships between the required satellite separations, 
the service-area geography, the protection ratios, and other system 
parameters. 

II. REQUIRED SATELLITE SEPARATIONS 

Consider the single-entry interference between two down-link 
satellite communications circuits. The up-link calculation has been 
shown to be dual, i.e., of precisely the same form [6], The geometry is 
shown in Figure 1. The following notation is used: S - satellite, 

E - earth station, W - wanted network, I - interfering network, T - 
transmit, R - receive. These symbols will also be used as subscripts in 
the equations below. It should be noted that the angle ^ is a 
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(a) Overall geometry 


SW 0 SI 



(b) Detail with Earth radius exaggerated for clarity 


Figure 1. Interference geometry between down-link networks. 
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two-dimensional vector since, for elliptical or shaped beams, not only 
its magnitude is important, but also the orientation of its plane with 
respect to the plane defined by the beam axis and the beam-maximum (or 
other reference) direction, e.g., the ellipse major axis for elliptical 
patterns. Similarly the angle ^ is a vector, but ^ can be treated as 
a scalar since there is no incentive for earth stations to use 
non-circular beams. 

The carrier and interference powers can be determined by means of 
the Friis transmission formula [7] and combined to give a 
well -approximated single-entry carri er-to-interference ratio 

E SWT °swt<V 6 swt> 

(C/Dcur = p + (1) 

SIT ° SIT ('l'2> G SIT^ n EWR^3 ,G EWR ) 

where E denotes effective isotropic radiated power, n antenna 
discrimination relative to the beam maximum, and G antenna gain in the 
beam-maximum direction [8], For satisfactory performance the 
carri er-to-interference ratio must equal or exceed the required 
protection ratio, which is the product of a co-channel protection ratio 
P and a relative protection ratio p(f), where f denotes the frequency 
offset from co-channel [9]. Therefore Equation (1) shows that the 
minimum allowable satellite spacing is implied in 

R DN = °SIT^2» G SIT^ °EWR^3’ G EWR^ ^ 
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where 


R DN P E SWT E SIT °SWT ^1* G SWT^ * 

The first four factors in are known system parameters. Also since 
calculations will always be performed at test points on the boundary of 
a service area and since, in practice, satellite beams will be shaped to 
give a*. reduction of approximately 3 dB at these test points, one can set 
d SWt( 1 1'1»^SWt) “-1/2. The left side of Equation (2) can therefore be 
considered known in an orbit synthesis procedure. 

It is important to note that Equation (2) is an implicit equation 
relating the required satellite separation to the separation of the two 
service areas. The existence of such a relationship has long been 
recognized qualitatively [10,11]. Equations (2) and (3) state 
quantitatively that, when frequency isolation is insufficient, the 
system requires antenna discrimination to achieve a required protection 
level, and that this is obtained as the product of the earth and 
satellite antenna discriminations. 

III. SEPARATIONS FOR CIRCULAR BEAMS 

For circular beams, the angle ^2 in Equation (2) becomes scalar and 
it is possible to solve explicitly for ^3 as a function of ^2 when the 
discrimination patterns Dsn> n EWR are specified. The relationship can 
be plotted conveniently as a universal set of contour curves, with Rq N 
as parameter and normalized values of 'i>2»’l>3 as coordinates. The need 
for normalization is implied in Equation (2) by the presence of the 
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gains G as arguments in the antenna discriminations. For example, when 
Gewr is large, corresponding to a narrow Earth -station beam, a given 
value of De^r is reached with a smaller value of ^3 than when Gewr is 
small. The universal curves are shown in Figure 2 for discrimination 
pattern envelopes recommended by the CCIR for co-pol ari zed FSS 
antennas [12,13]. Unfortunately, no corresponding cross-polarized 
patterns have been recommended as yet. Two sets of curves are required 
because of the piecewise specifications of 0 (^ 3 ). At first sight it 
might appear that four sets would be required since four "pieces" are 
used in the CCIR specification of Earth-station discrimination. 

However, it turns out that the far side-lobe and back-lobe regions imply 
so much discrimination that the interference can be neglected, while the 
constant part of the pattern, representing the near-sidelobe envelope, 
is accounted for by the discontinuity in the ^3 values which results 
from switching between Figures 2a and 2b in accordance with the 

directions given below them. 

The expression of the universal curves in terms of the 

antenna-centered "off-axis" angles ^ and 4*3 is natural and also useful: 
for example, it shows that the "separation" of service areas, measured 
by i|» 2 * is rigorously the distance from the aimpoint of the interfering 
satellite antenna to that test point of the area suffering interference 
which is on the highest discrimination contour of the interfering 
satellite antenna pattern on the Earth surface. For circular beams and 
high satellite elevations this would be the test point nearest to the 
interfering satellite aim point. Nevertheless, for many system 
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specified G^ WR means the numerical gain value. 



Figure 2 . Universal curves for the minimum al low-able satellite spacing 
angle 1P3 as function of the normalized off-axis angle ^ 

<p is the half-power beam width of the satellite antenna) d/A 
tne diameter-to-wavelength ratio of the EWR antenna. 
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calculations it is more useful to find the geocentric satellite 
separation A<|> instead of the topocentric angle directly as a 
function of the longitude differences and latitudes of EWR, EIR, and SI 
instead of ^ Actually the angle ^ turns out to be a rather good 
approximation for A<J>, which can be improved sufficiently for all 
practical purposes by 

A 4> = <|> 3 [1.023-0.302 cos(<j> M -<f> E )cos0 E ] 1/2 rad (4) 

where is the longitude of the midpoint between the wanted and 
interfering satellites (SW,SI) and 4> E , 0 E are the longitude and latitude 
of the test point (EWR), respectively. 

The relationship between i|> ? and the geocentric variables can be found 
by substituting suitable expressions obtained from Figure 1 into the 
cosine law 

cosily = (b^+c^-g^)/2bc , (5) 

but for many purposes the relationships 
9 9 

6.54R e < *2 < 5.62R e rad * (6) 

where R^ is the earth radius, give a sufficiently good estimate. The 
correct value is close to the upper limit for high satellite SI 
elevations and to the lower for low elevations, as viewed from the aim 
point EIR. The variation of the required separation A<f> for various 
system parameters and configurations in terms of longitude and latitude 
is shown in Figures 3 and 4. From these and more such computations the 
following results emerge [14]: 


7 



NIN. GEOCENTRIC SPRCING IOEG.I 


ORIGINAL PAGE IS 

OF POOR QUALITY 



(a) IWP» AT 0*N 



JM EWR AT E0‘ N 


(a) 



(cl EWR AT 40* N 


(b) 



(c) 


(d) 


Figure 3. Minimum geocentric satellite spacing when earth stations are 
separated in longitudinal direction. = 35 dB, G$ij = 40 
dB, G^WR = 50 dB. 
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(a) 


(b) 



Figure 4. Minimum geocentric satellite spacing when earth stations are 
separated in latitudinal direction. Rq^j = 35 dB, 6517 = 40 
dB , SgviR = 50 dB . 
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(a) the smallest required separation occurs for practical geometries 
when the wanted satellite is near the longitude of the center of 
its service area, 

(b) for a substantial range of orbital locations about this longitude 
the required separation varies little. 

This last result, which appears to be true also for elliptical beams 
(see Section IV), is very important in the synthesis procedure discussed 
in Section V, because it reduces or eliminates the need to recalculate 
the required satellite separations as satellite orbit assignments are 
changed. 

IV. SEPARATIONS FOR ELLIPTICAL BEAMS 

For elliptical beams, the required satellite separations can be 
calculated by the following procedure: 

(1) For each service area choose test points on the periphery and 
calculate the orientation and axial ratio A r of the minimum 
projected ellipses which enclose them, using 2° increments in 
satellite longitude [15], 

(2) Choose a wanted and an interfering satellite and a longitude 
within the orbital arc under consideration, and consider the 
desired and interfering satellites collocated at that longitude. 

By incrementally increasing A<(>, with the satellites located 
symmetrically on each side of the initially chosen central 
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longitude, calculate the required satellite separation from 
Equation (2) with 

% 

•¥ ^ 9 9 O 1 /O 

D SIT ( ^2» G SIT ) = D s I T C 7" (C0S 9 + A r Sin 9)1 » (7) 

where i}> 0 is the half-power beamwidth in the direction of the major 
axis and 6 is the angle between the plane in which ^ is measured 
and the plane determined by the SIT beam axis and the SIT ellipse 
major axis. Equation (5) is used for evaluating ^ in this 
calculation. 

(3) Repeat step (2) for other central longitudes. The result for the 
worst test point for each service area shown in Figure 5 is shown 
in Table 1 for a sequence of central longitudes. For each pair of 
satellites two separation values can be obtained, depending on 
which service area is considered protected, with the satellite of 
the other area considered interfering. The larger of these 
appears as A<j> in Table 1. Notice the slow variation of A<|> with 
the mean satellite longitude, provided that this longitude is not 
too different from that of the protected service area. One 
satellite per service area was assumed in these calculations. 

(4) Choose the largest of all the A<J> values obtained in this process 
over the allowable orbital arc for each satellite pair and denote 
it As. The triangular matrix of As elements corresponding to 
Table 1 is shown in Table 2 for all satellites constrained to the 
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SERVICE AREA 

ARGENTINA 

BOLIVIA 

CHILE 

PARAGUAY 

PERU 

URUGUAY 



Figure 5. Geography of the six-service-area scenario. Dots indicate 
test points. 
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TABLE 1 


A<j> VALUES IN DEGREES 

Single-entry protection ratio = 30 dB. The discrimination 
patterns of refs. 16,17 were used for the Earth antennas with 
d/X = 60 and 0.6 aperture efficiency. Those of ref. 18 were 
used for the satellite antennas with beam width as determined 
by the minimum ellipse. It was assumed that SW and SI produce 
the same power density at their respective aim points, 
corresponding to « 33 dB. 


Satell ite 
pai rs 

Mean Satellite Longitude 

70°W 

80°W 

90°W 

100°W 

110°W 

1 

ARG 

BOL 

4.00 

4.02 

4.05 

4.12 

4.17 

ARG 

CHL 

4.18 

4.05 

4.00 

4.02 

4.19 

ARG ! 

PRG 

4.24 

4.28 

4.32 

4.28 

4.32 

ARG ' 

PRU 

0.94 

1.04 

1.15 

1.25 

1.41 

ARG 

1JRG 

4.18 

4.14 

4.06 

4.06 

3.94 

BOL 1 

CHL 

4.13 

4.20 

4.28 

4.39 

4.57 

BOL 

PRG 

4.00 

3.99 

4.00 

4.03 

4.04 

BOL ' 

PRU 

3.87 

3.95 

3.99 

4.10 

4.26 

BOL ' 

URG 

0.39 

0.38 

0.38 

0.84 

0.94 

CHL 

PRG 

1.08 

1.14 

1.25 

1.46 

2.00 

CHL ' 

PRU 

3.84 ' 

3.83 

3.85 

3.89 

3.94 

CHL ' 

URG 

0.42 

0.41 

1.03 

1.28 

1.59 

PRG 

PRU 

0.49 

0.49 

0.50 

0.85 

1.10 

PRG 

URG 

2.16 

2.19 

2.20 

2.34 

2.46 

PRU | 

URG 

0.40 

0.37 

0.37 

0.33 

0.32 


13 


c> 


80°W to 110°W orbital arc. These separations will be used as 
constraints in the synthesis procedures of the next section. 


TABLE 2 

AS VALUES IN DEGREES 



BOL 

CHL 

PRG 

PRU 

URG 

ARG 

4.17 

4.19 

4.32 

1.41 

4.14 

BOL 


4.57 

4.04 

4.26 

0.94 

CHL 



2.00 

3.94 

1.59 

PRG 




1.10 

2.46 

PRU 





0.37 


V. LINEAR PROGRAMMING FORMULATION 

The FSS synthesis problem can now be formulated as a linear program 
with a set of nonlinear side constraints. The set of satellite 
locations which satisfy the constraints constitutes the feasible region. 
A variety of linear functions can be selected to be optimized. Three 
functions have occurred to us: 


(a) to search only for some point in the feasible region by setting 
the function to be minimized equal to zero; 

(b) to minimize the occupied orbital arc; 

(c) to minimize the sum of the absolute deviations of the satellite 
locations from a specified set of locations. 


14 



The last objective has been implemented in the formulation described 
below. The following parameters and non-negative variables are used: 

ej(wj) - easternmost (westernmost) feasible location for satellite 
j, in °W, 

dj - desired location for satellite j, 

As-jj - required minimum separation between satellites i and j, as 
in Table 2, 

xj - relative location (in degrees west of ej) of satellite j, 

Xj + (xj-) - degrees west (east) of dj that satellite j is located, 

P i j ( n i j ) “ de 9 rees west (east) of satellite j that satellite i is 
located (defined for i < j). 

The FSS orbit assignment synthesis is then formulated as follows: 


Minimize 

2=1 (Xj + + Xj~) (8) 

subject to 

Xj " x j + + x j“ = d j for all j , (9) 

Xj - Xj - pjj + n-j j = 0 for all i,j 

where i < j , (10) 

Pij + n ij > ASi j for all i,j 

where i < j , (11) 
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x j < w j - e j 

for al 1 j , 

(12) 

X J» x J + » x j" > 0 

for all j , 

(13) 

Pi j * n ij > 0 

for all i , j 



where i < j , 

(14) 

Pij ' "ij = 0 

for all i , j 



where i < j . 

(15) 


The objective function. Equation (8), computes the sum of the 
absolute deviations of the prescribed satellite locations (xj's) from 
the desired or ideal locations (dj*s). These absolute deviations are 
measured in the first set of constraints. Equation (9). The actual 
separations between all pairs of satellites (pij + n-jj) are computed in 
the second set of constraints. Equation (10), and are compared to the 
minimum required satellite separations in the third set of constraints. 
Equation (11). The constraints of Equation (12) guarantee that the 
location prescribed for each satellite is feasible. Constraint Equations 
(13) and (14) indicate that all of the variables in the problem are 
restricted to non-negative values. Finally the complementary 
relationships between pairs of variables p-j j and n-jj are enforced by the 
constraints of Equation (15). 

Linear programs are much more readily solvable than nonlinear 
programs and integer programs. They are most often solved by the simplex 
method [19], This technique examines a sequence of basic solutions to 
the constraints of the linear program. Each solution examined has an 
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objective function value no less favorable than that of the previous 
solution. The algorithm terminates when it is determined that no 
improved solution can be found. 

The presence of the nonlinear side constraints of Equation (15) 
prevents us from using the simplex method in its most common form. The 
method can be modified to handle these additional constraints through the 
use of restricted basis entry: p-jj can not be a basic variable if n-fj is 
a basic variable, and vice versa [20]. When employing the simplex method 
with restricted basis entry, we are certain to find a local, but not 
necessarily a global, optimum. As formulated, the problem has m(m+2) 
variables, where m is the number of satellites, and m^ constraints, not 
counting the simple bound constraints of Equations (12) to (14) and 
complementarity constraints of Equation (15). The formulation is similar 
to one suggested by Ignizio for the N-job, single-machine scheduling 
problem [21], 

VI. MIXED INTEGER PROGRAMMING FORMULATION 

The same problem can also be formulated as a mixed integer program 
[22], A global optimum is guaranteed when this formulation is employed. 
However, the computational effort required to find a final solution can 
be immensely greater than it would be with the linear programming 
formulation. In any case, this formulation is helpful in assessing the 
quality of the solutions found using the linear programming formulation 
on small test problems. 
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To complete this formulation, we need the following definitions: 


' = min {e,} 

J J 

9 

(16) 

- = max {w.} 
j 4 

9 

(17) 

1 if 

satellite i is located west of 


..j = { satellite j. 

(18) 

0 otherwise . 



The objective function and the constraints (9), (10), (11), (13) and 
(14) appear in this formulation precisely as they did in the linear 
programming formulation. The complementarity constraints (15) are 
replaced by two new constraints: 

(e-w) x i j + p.j < 0 (15a) 

(w-e)x.. j + n i j < (w-e) (15b) 

which, together with the nonnegativity restrictions on pij and n-j j , 
constraints (14), guarantee that either pjj = 0 or n^ j = 0 for all pairs 
of satellites i and j. 

If there are m satellites, the mixed integer formulation entails 
m(m+2) continuous variables, m(m+l)/2 binary variables, and 2m 2 -m 
constraints specified by Equations (9), (10), (11), (15a), (15b). The 
time required to solve an FSS synthesis problem with this formulation 
will be most heavily dependent upon the number of binary variables. For 
large problems (many satellites), this formulation may involve 
prohibitive solution times. 
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VII. NUMERICAL RESULTS 

The FSS synthesis minimizing the sum of absolute deviations of 
orbital positions from a prescribed "desired" set was solved, both as a 
linear program with the simplex method with restricted basis entry and 
as a mixed integer program via branch-and-bound [23]. The service areas 
and test points were those of Figure 5 with one satellite per service 
area. The available orbital arc for each satellite was specified as 
80°W to 110°W. It was assumed that each satellite would carry a full 
complement of frequency channels, so that a co-channel calculation is 
appropriate. A single-entry C/I value of 30 dB was chosen with the 
intent of achieving a 25 dB aggregate co-channel C/I ratio. With these 
assumptions the As values of Table 2 are pertinent. Three problems were 
run, differing only in the specified "desired" satellite locations. In 
problem 1, this "desired" location was specified for every satellite as 
95°W, the center of the arc. In problem 2, all "desired" locations were 
specified at 110°W, the westernmost end of the arc. In problem 3, each 
was specified near the central longitude of the ellipse circumscribing 
the service area to be served; these "desired" longitudes are indicated 
in the column labeled DL in Table 3, which shows the solutions obtained 
for all three problems by both methods. The LP formulations required 48 
variables and 36 constraints, while 63 variables, 15 of them binary, and 
66 constraints were needed for the MIP formulation. 

The solutions to these test problems illustrate some important 
points. First of all, the solution of a synthesis problem by means of 
an integer program can require a substantially greater amount of 
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computer time than by means of a linear program. Secondly, the two 
approaches used can produce strikingly different solutions. (See the 
results for Problems 1 and 2.) It may also happen that the same 
solution will be found with both methods, even though this is not 
evident from the results presented here. Finally, acceptable solutions, 
in terms of aggregate co-channel C/I ratios, are obtained even when the 
objective function value for the linear programming solution differs 
substantially from the mixed integer programming solution, the global 
optimum. This is not unexpected because the as constraints guarantee 
acceptable single-entry C/I ratios. Table 4 shows the distributions of 
aggregate co-channel C/I ratios for the two methods and three problems. 
It will be remembered that a 30 dB single-entry constraint was used to 
calculate the As table on which all these calculations are based. 

TABLE 4 

NUMBER OF TEST POINTS CORRESPONDING 
TO A GIVEN AGGREGATE CO-CHANNEL 
C/I RATIO RANGE FOR EACH PROBLEM AND METHOD 


Problem 

Method 


<27 

C/1 

27-28 

Interval 

28-30 

wmm 

EUB&l 

>35 

1 

LP 


0 

1 

8 

16 

29 

1 

MIP 


0 

1 

6 

18 

29 

2 

LP 


0 

0 

5 

20 

29 

2 

MIP 


0 

4 

10 

25 

15 

3 

LP 


0 

0 

4 

16 

34 

3 

MIP 


0 

0 

9 

14 

31 
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VIII. CONCLUSIONS 


An implicit relationship has been derived which relates the 
topocentric separation of two satellites as required for a given level 
of single-entry protection to the separation and orientation of their 
service areas. For circular beams and topocentric angles the results 
are presented explicitly; a computational approach is given for 
elliptical beams and for use with longitude and latitude variables. It 
is found that the geocentric separation depends primarily on the service 
area separation, secondly on a parameter which characterizes the 
electrical design, and only slightly on the mean orbital position of the 
satellites. Both linear programming (LP) and mixed integer programming 
(MIP) algorithms have been implemented, with the sum of the absolute 
deviations of the orbital locations from a prescribed "ideal" set as 
objective function. Three "ideal" sets were used. A single-entry 
protection ratio of 30 dB was specified with the intent of satisfying an 
aggregate co-channel C/I ratio of 25 dB. The set of orbital locations 
chosen by the LP and MIP methods differed substantially in two cases 
and, to a lesser degree, in the third. Still, all six solutions which 
were found provided acceptable protection ratios at all 54 test points. 
The worst aggregate co-channel C/I ratio found was 27 dB. 

The results are encouraging with respect to applying the LP 
procedure to larger scenarios. The MIP formulation may result in 
excessive computation times when many satellites are involved, but it is 
guaranteed to arrive at a global optimum; it will therefore be useful 
for evaluating the efficacy of the LP approach via smaller test 
problems, such as the three presented here. 
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